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Abstract 

The  optical  field  in  free  electron  lasers  can  sometimes  develop  side- 
bands and  exhibit  very  complex  behavior.  The  purpose  of  this  report 
is  to  develop  a  two  mode  model  of  the  optical  field  in  a  free  electron 
laser  in  which  the  optical  modes  are  coupled  by  the  free  electrons  in 
the  laser.  This  model  may  be  used  to  study  the  long  term  stability 
of  the  optical  wave  and  the  development  of  sidebands  in  the  optical 
field.  In  this  report  the  equations  of  motion  governing  the  interaction 
of  an  electron  with  the  optical  modes  will  be  derived. 

1.  Introduction 

Initial  theoretical  studies  of  the  free  electron  laser  (FEL)  used  quantum  me- 
chanics ([l]-[3]).  However  it  was  shown  in  [4]  that  a  classical  treatment  of 
FELs  would  adequately  describe  most  phases  of  operation  ([5]-[10])  and  re- 
sult in  substantial  simplification  of  the  resulting  equations  describing  their 
operation. 

One  way  to  understand  some  of  the  main  features  of  the  FEL  problem 
is  to  consider  the  motion  of  a  pendulum  with  weak  damping  and  a  nearly 
constant  tangential  forcing.  In  a  mathematical  model  of  the  FEL  with  a 
tapered  magnetic  wiggler,  however,  the  weak  damping  term  in  the  pendulum 
equation  cannot  be  ignored  because  the  wiggler  is  extremely  long  compared  to 
its  period,  and  one  must  account  for  the  "cumulative"  effects  of  this  damping 
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over  the  entire  length  of  the  wiggler  in  order  to  obtain  accurate  solutions. 
The  only  way  to  systematically  account  for  these  higher  order  effects  is  to 
formulate  a  completely  dimensionless  problem  which  contains  dimensionless 
parameters  whose  magnitudes  characterize  various  distinguished  limits. 

In  a  FEL  frequently  the  optical  field  can  develop  sidebands  and  become 
highly  complex  in  its  behavior.  The  purpose  of  this  report  is  to  develop  a 
model  in  which  there  are  only  two  modes  present  in  the  optical  field:  these 
modes  are  coupled  to  one  another  through  the  electrons  in  the  FEL.  This 
model  can  be  used  to  study  the  long  term  stability  of  the  optical  wave  and  the 
development  of  sidebands  in  the  optical  field.  A  regime  of  particular  interest 
is  that  of  sustained  resonance  in  the  two  mode  model.  In  this  report  we  shall 
derive  the  equations  of  motion  governing  the  interaction  of  an  electron  with 
the  optical  fields  in  the  two  mode  model  of  a  FEL. 

2.1  Equations  of  Motion 

We  shall  formulate  the  equations  of  motion  as  a  first  order  system  of 
ordinary  differential  equations.  Our  derivation  follows  that  in  [11].  However 
we  shall  be  extending  the  derivation  given  there  for  a  single  optical  mode 
to  the  two  optical  mode  case.  The  convention  we  shall  use  is  that  upper- 
case letters  will  represent  dimensional  variables  and  lower-case  letters  will 
represent  dimensionless  variables. 

To  determine  the  leading  order  equations,  consider  a  single  electron  mov- 
ing in  a  helical  wiggler  magnetic  field  given  by 

Bw  =  Bw    X  cos  —  Z  +  Y  sin  —  Z  L  (1) 

I  Aw  Aw       ) 

where  Bw  is  the  magnetic  field  amplitude  and  \w  is  its  wavelength.  The 
vectors  X,  Y,and  Z  are  orthonormal  and  are  associated  with  the  coordinates 
X,  Y,  and  Z,  where  Z  is  the  coordinate  along  the  direction  of  the  magnet  axis. 
We  shall  assume  that  the  radiation  (or  signal)  field  is  a  linear  superposition 
of  two  modes  and  is  given  by 

E»  =  Ei1)  +  E<2\  (2) 

where 

E^  =  -^(Z){Xsin0t  +  Ycos0i},i  =  1,2  (3) 


and 

Ba  =  Bi1)  +  B?>,  (4) 

where 

B»  =  ZxE«  t  =  l,2.  (5) 

E^  is  the  radiation  amplitude  of  mode  i  and 

0t  =  ^Z-^')r  +  ^),z  =  l,2  (6) 

As 

where  A^  is  the  wavelength,  u^  is  the  frequency,  T  is  time,  and  <fty  is  the 
phase  of  mode  i,  i  =  1,2. 

We  shall  assume  that  the  electron  is  relativistic,  though  its  transverse 
velocity  (the  velocity  orthogonal  to  the  Z  axis)  is  small.  Let  (3  be  the  dimen- 
sionless  electron  velocity  (/3  is  the  electron  velocity  divided  by  c,  the  speed 
of  light).  We  shall  decompose  /3  into  a  component  parallel  and  a  component 
perpendicular  to  the  Z  axis: 

/9   =  fl,Z  +  /J±,       0< /9±  «  0j|  <  1.  (7) 

The  equation  of  motion  (Newton's  second  law)  for  the  electron  is 

P  =  eEs+e/3  x(B5+Bu;)+space  charge  force-f  radiation  reaction  force,  (8) 

where  e  is  the  electron  charge  and  P  =  7mc/3    is  the  electron  momentum.  7 
is  the  Lorentz  factor,  defined  by 

y-?TTT'  (9) 

c  is  the  speed  of  light,  and  m  is  the  electron  rest  mass.  Criteria  for  neglecting 
the  effects  of  the  radiation  reaction  force  and  the  space  charge  force  were  given 
in  [12],  and  were  derived  by  examining  higher  order  terms  in  the  equations 
of  motion.  These  conditions  will  be  given  later  in  this  report. 
The  conservation  of  energy  implies  (7771c2)  =  ec/3  •  Es  or 

7  =  —  fi     V*  (10) 

mc 

To  gain  a  rough  idea  of  the  approach  to  be  followed,  we  shall  assume 
initially  that  only  a  single  mode  in  the  radiation  field  is  present.    Later  we 


shall  repeat  the  same  calculations  in  more  detail  with  two  modes  in  the 
radiation  field. 

The  momentum  equation  (8)  implies 

(7mc/3  )  =  7(mc/3  )  +  7mc/3   =  73(mc/3  )(/3  •  /3  )  +  -ymc/3  ,  (11) 

since  7  =  73/3   •  /3  . 

Combining  the  momentum  equation  (11)  and  the  energy  equation  (10), 
we  conclude  that 

7mc/3   =  eEs-e/3(/3  •  E,)  +  e0   x  (Z  x  E,)  +  e(3   xBw.        (12) 

Since  /?   x(ZxEs)  =  Z(/3  •  E,)  -  E,(/3  •  Z),  equation  (12)  may  be  written 
as 

7mc/3   =e(l-/?n)E3  +  eZ/3   •  Es  +  e/3    xB„  (13) 

Using  the  fact  that  1  —  /3u  <C  1  and  that 

0  -E8  =  /3x.Ea   ,#j.  <1,  (14) 

we  obtain  the  approximate  equation 

(3  =  — 0   x  B„.  (15) 

7mc 

Substituting  /3  ,  the  solution  of  equation  (15),  into  the  energy  equation  (10), 
we  find 

7  = sin  0,  (Id) 

mc 


2ttc  /        Xw  1  +  bj 
K   V         K     27 
where  <f>  =  (2k /\w  +  2n/\3)Z  —  ujsT  +  <^>5,  and 


*  =  ^F  ( 1  -  T^1^    .  (17) 


4.  =  5^,  (18) 

is  defined  as  the  dimensionless  amplitude  of  the  wiggler  in  Section  (2.2).  The 
details  of  the  above  calculations  (for  the  case  of  two  modes)  will  be  given  in 
Section  (2.3). 


We  can  now  define  a  (dimensionless)  resonant  electron  energy  7,.  for  which 
(j>  =  0.  This  energy  is  obtained  from  equation  (17)  and  has  the  value 

^  ,*-!+&.  (19) 


A 


A  free  electron  laser  usually  operates  in  a  state  for  which  7  «  7,.  (see  [2],  [6], 
and  [7]).  Equation  (16)  indicates  that  electrons  with  (dimensionless)  energy 
7r  in  a  magnetic  wiggler  with  wavelength  A^  and  field  amplitude  bw  will 
spontaneously  emit  radiation  (if  7V  <  0)  or  absorb  radiation  (if  7V  >  0)  of 
wavelength 

A.  =  Ml+£).  (20) 

2.2  Normalization  and  Nondimensionalization 

It  is  very  important  to  decide  on  characteristic  scales  for  the  dependent 
and  independent  variables  in  the  study  of  a  system  of  differential  equations 
by  asymptotic  methods.  The  appropriate  choice  of  scales  is  usually  dictated 
by  the  physical  regime  of  interest  represented  in  the  system  under  study  and 
can,  in  fact,  be  changed  to  suit  a  different  regime.  They  key  point  is  that 
once  a  set  of  scales  is  chosen  all  dependent  and  independent  variables  become 
dimensionless  and  order  one  in  magnitude.  There  then  is  no  ambiguity  about 
the  size  or  numerical  magnitude  of  the  dependent  and  independent  variables. 
Of  course,  there  are  now  dimensionless  parameters  present  in  the  problem 
which  are  functions  of  the  characteristic  scales.  It  is  their  orders  of  magnitude 
which  indicate  what  physical  regime  is  under  study,  and  what  is  small,  etc. 

We  shall  use  Xw/2tt  =  L0  as  a  characteristic  unit  of  length  and  L0/c  =  T0 
as  a  characteristic  unit  of  time.  Let  27r/Aw  =  Kw  be  the  wavenumber  of  the 
wiggler.  Then  in  terms  of  these  units  we  define  the  dimensionless  variables 

z  =  KwZ  =  ^-   ,    t  =  KwcT  =  ^.  (21) 

(Note  that  lower  case  variables  will  always  be  dimensionless.) 
We  now  introduce  the  dimensionless  field  variables: 

bw  =  1r^—1Bw,    fcf  =  — L-jfl„    ea  =  -±—Es,  (22) 

Kwmci  Kwmcl  Kwmcz 


where  mc2  =  0.511  MeV.  Using  the  normalizations  introduced  in  equations 
(21)  and  (22),  equations  (8)  and  (10)  may  be  written  in  dimensionless  form 

d(l(3) 


dt 


=  ea  +  /?    x  (bs  +  bw)  +  higher  order  terms,  (23) 

J  =  0-e..  (24) 

In  equation  (23),  the  higher  order  terms  include  the  normalized  space  charge 
force  and  the  normalized  radiation  force. 

To  characterize  the  order  of  magnitude  of  each  term  in  these  equations 
we  need  to  introduce  a  small  parameter  £,0  <  £  <  1.  One  way  to  do  this 
is  to  choose  the  parameter  £  to  be  a  measure  of  the  slow  variation  of  the 
wiggler  field  along  the  magnet  axis.  To  do  this  we  define  e  as 

*w         Lq  1 

£  =  2^I  =  T'  or  £  =  1  (25) 

where  L  and  1  are  respectively  the  dimensional  and  dimensionless  length  of  the 
wiggler  and  L0  is  the  characteristic  length  defined  by  the  wiggler  in  equation 
(21).  The  smallness  of  £  signifies  that  the  magnetic  wiggler  produces  a  small 
effect  which,  for  an  electron  moving  along  the  magnet  axis,  lasts  for  a  long 
time. 

Using  the  small  parameter  e  we  can  define  a  "slow"  variable  z  as 

z  =  ez,     0<£<1.  (26) 

The  slow  variable  z  enables  us  to  define  "slowly  varying"  functions  of  z  in  the 
following  sense:  an  0(1)  function  of  z  which  varies  slowly  along  the  wiggler 
axis  while  maintaining  its  order  of  magnitude  can  be  written  as  f(£).  Note 
that  df/dz  =  0(e).  A  slowly  varying  function  f  of  the  dimensionless  variable 
z  may  be  written  as 

/(i)  =  /(«)  =  /(^£)  =  /(f), 

so  that  slowly  varying  means  changes  occur  over  length  scales  which  are  much 
smaller  than  the  length  of  the  wiggler. 

Typically  the  wiggler  length  L  is  several  meters  and  the  wavelength  Xw 
of  the  wiggler's  magnetic  field  is  on  the  order  of  a  centimeter,  so  e  is  a  small 


number,  about  10~3.  Also  since  the  dimensionless  wiggler  length  1  is  1/e,  the 
wiggler  must  be  considered  long  in  the  mathematical  sense.  This  means  that 
the  cumulative  effects  of  higher-order  terms  may  be  important. 

We  will  now  denote  all  normalized  (dimensionless)  variables  with  an  over- 
bar.  The  wiggler  magnetic  field  is  then  written  in  normalized  form  as 

bw  =  bw(z){Sicosz  +  ysin2}.  (27) 

It  can  be  seen  now  from  equation  (19)  that  bw  ^>  1  leads  to  7;?  ^>  ^-,  ie., 
extra  high  energy  electrons  are  required.  On  the  other  hand,  if  bw  <C  1, 
the  wiggler  field  is  then  too  weak  to  significantly  decelerate  electrons.  We 
conclude  that  the  only  reasonable  choice  is  bw  =  bw(z)  =  0(1).  Because  the 
operation  of  free  electron  lasers  is  near  the  resonant  state  defined  in  equation 
(19)  (see  references  [2],  [6]  ,  [7]),  and  since  typical  values  of  7r  are  large, 
we  shall  write  the  dimensionless  energy  of  an  electron  and  the  dimensionless 
wavenumber  of  the  signal  in  the  following  form  consistent  with  equation  (19): 

7  =  £"o/27,  (28) 


-f  =  z~°~k-  (29) 

In  the  above  equation,  ka  is  the  0(1)  dimensionless  signal  wavenumber 
and  ea  =  n  is  a  second  small  parameter  which  measures  the  ratio  \s/\w. 
The  determination  of  /z  comes  from  a  consideration  of  the  physical  regime 
under  discussion.  Recall  that  variables  with  an  overbar  are  dimensionless 
and  0(1). 

After  this  preliminary  analysis,  let  us  assume  again  that  there  are  two 
modes  present  in  the  radiation  field.  In  normalized  form  its  components 
appear  as 

eW  =  -£*§(')  =  {xsin  0{  +  y  cos  0,},  (30) 

b<*>  =  z  x  e<*>,  (31) 

for  i=  1,2,  and  where 

Si  =  e-"V;\z  - 1)  +  J®,  (32) 

again  for  i  =  1,2.   Note  that  the  total  radiation  field  is  given,  as  before,  by 
equations  (2)  and  (4).  The  parameter  (  in  equation  (30)  is  again  determined 


by  the  physical  regime  under  discussion.  We  shall  show  in  section  [2.3]  that 
the  only  really  interesting  limit  for  the  free  electron  laser  is  £  =  0. 

The  radiation  parameters  ey'  and  <j>y  depend  on  the  wave  evolution  pro- 
cess and  must  be  determined  by  Maxwell's  equations.  There  are  some  open 
questions  here  which  will  be  discussed  later.  We  will  see  that  es  evolves  on 
a  slower  scale  than  bw. 

There  are  basically  three  types  of  magnetic  wigglers  or  undulators.  There 
are  those  in  which  A^,  varies,  those  in  which  bw  varies,  and  those  in  which 
both  Au,  and  bw  vary.  Our  approach  to  this  problem  can  be  adapted  to  any 
of  these  wigglers.  However,  the  work  in  [7]  suggests  that  the  constant  \w 
case  is  probably  the  best  method  for  enhancing  efficiency.  Since  it  also  offers 
apparent  hardware  advantages  we  will  restrict  our  attention  to  this  case. 

2.3  Higher  Order  Correction  Terms 

In  this  section  we  calculate  higher  order  correction  terms  for  the  equations 
of  motion  and  we  mention  conditions  derived  in  [12]  for  neglecting  the  space- 
charge  force  and  the  radiation  reaction  force. 

We  begin  by  looking  at  equation  governing  the  temporal  evolution  of  the 
perpendicular  component  of  the  velocity,  0±.  From  equations  (23)  and  (24) 
we  find 

e-a/2jt(P±)  =  es  +  (/?„*  +  0±)  x  (z  x  e8  +  bw).  (33) 

Since  e8  is  perpendicular  to  z,  we  find  that  the  last  equation  may  be 
written 


t-a,2  —  {^)  =  (1  -  /?||)es  +  /?||5  x  bw  +  (  higher  order  terms)  (34) 

Now  ,  for  a  single  electron,  4f  =  /Sjj.  We  can  regard  this  as  an  ordinary 
differential  equation  expressing  z  as  a  function  of  t.  A  little  later  it  will  be 
convenient  to  invert  this  relationship  and  to  express  t  as  a  function  of  z. 

From  [12],  the  following  conditions  must  hold  in  order  to  neglect  the 
radiation  force: 

(^)  =  o(e-<)  (35) 

Conditions  necesssary  for  neglecting  the  space  charge  force  are 


^?L 


=  o(cba'^)  (36) 


Once  a  is  determined  (and  it  will  be  shown  later  that  £  is  zero),  the  above 
two  equations  provide  conditions  which  must  hold  in  order  to  consistently 
neglect  the  effects  of  the  radiation  and  space  charge  forces  on  an  electron.  In 

the  second  equation,  the  ratio  -§•  is  the  ratio  of  the  plasma  frequency  to  the 
signal  frequency.  It  will  turn  out  that  this  ratio  is  a  fundamental  parameter 
measuring  the  rate  of  evolution  of  the  radiation  field.  Consequently  the 
second  condition  on  the  neglect  of  the  radiation  force  plays  an  important 
role  in  the  evolution  of  the  optical  field. 

Inserting  the  specific  field  quantities  into  equation  (34)  then  gives 

e_0f/2-j7(/?±)    =    -P\\bw(z)  (x  sin  z  -  y  cos  z) 

-(1  -  P\\)&M  (x  sin  0!  +  y  cos  6X)  (37) 

-(1  -  /?||)e*&<2)  {x  sin 02  +  y  cos 02) 
-f  higher  order  terms, 

where  b\  =  e\  since  b| ,    =  z  X  el  . 

Noting  that  0u  dt  =  dz,  we  integrate  equation  (37)  to  obtain  the  form 

/?x     =    e      ■—  (£  cos  z  4-  y  sin  z) 
1 
v 

_e3a/2_j£  /£  gjn  z  _   *  CQS  z\ 

1 

-c3°/2+(iT^r)  (*  cos  9i -$  sin  *i)      •  (38) 

1(2) 

-e3o/2+^r%  (* cos  h-y sin  fc) 

+0(e5o,/2) 

where  0,  =  e~afc^(2  —  <)  +  <f>[^-  In  equation  (38),  b'w  stands  for  the  derivative 
of  bw  with  respect  to  eaz.  However  since  bw  is  a  function  of  z  =  e^r,  it  follows 


that 


db, 


=  0{tl-°) 


(39) 


d{e°z) 

The  0(e5Q/2)  error  in  equation  (38)  is  due  to  ignoring  the  variations  of 
dbw/d(tQZ),bs,  and  <f>B  with  z.   To  obtain  (38)  ,  we  have  also  assumed  that 
/?±(0~)  =  0(e5or/2),  ie.,  the  transverse  velocity  of  the  incoming  beam  can  be 
ignored  before  it  enters  the  magnetic  wiggler. 
Now 

(40) 


J  =  ^/3.es  =  ^.(e<1»  +  el2') 


since  the  component  of  the  electron  velocity  parallel  to  the  wiggler  axis  is 
orthogonal  to  the  signal  fields.  Inserting  equations  (37)  and  (38)  into  the 
above  expression  in  (40)  ,  we  find  after  some  calculation  and  rearrangement 
that 


drf 
~dt 


a+t 


T    T(l)  7    T(2) 

w-^~  sm{z  +  0T)  +  -^-  sin(z  +  T2) 


1 


1 


—  e 


2a+i 


VtO)  T'T(2) 

-^-  cos{z  +  0X)  +  -2J-  cos(z  +  02) 

7  7 


(41) 


+e 


2a+2^ 


i1.**.21 


gft?1 


-U=sin(#2-e,)  +  -^=sin(fl1  -«2) 
L7*<"  7tP 

+0(e3a). 

We  shall  now  define  the  phases  ^,-,i  =  1,2  by 

A :  =  z  +  Wi  =  z  +  e— Ii°(*  -  t)  +  #. 


(42) 


We  observe  that  since  many  of  the  parameters  in  the  problem  depend  on 
z,  as  previously  mentioned,  it  will  be  convenient  to  use  d/dz  rather  than  d/dt 
for  the  equations  of  motion.  The  velocity  dz/dt  can  be  determined  from  the 
relationships  dz/dt  =  0\\  and 


fl|  =  (l-l/72-#) 


1/2 


(43) 


10 


Using  the  above  equation  and  the  expressions  for  7  and  @±  in  equations 
(9)  and  (38)  ,  we  find  after  some  effort,  that 


■    .  1+,(!±S)t«..w(l±|l 


— e 


2a+Z 


(44) 


+0(eJa) 


Since 


it  follows  that 


d  7       cte  <f  7 
cfa         eft  c/z 


(45) 
(46) 


dj        1  c/7 
dz       /?||  eft 

We  thus  can  combine  equations  (41)  and  (44)  above  to  derive  an  expression 
for  the  rate  of  change  of  7  with  respect  to  the  dimensionless  distance  along 
the  wiggler  axis,  z.  We  find  after  some  calculation  that 


c?7 
dz 
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Equation  (47)  determines  the  evolution  of  7  with  respect  to  z,  and  it  can 
be  seen  that  this  depends  on  the  field  quantities,  and  on  the  phases  c/>i,c/>2 
defined  in  (42).  We  now  determine  the  evolution  of  these  latter  quantities 
with  respect  to  z.  Differentiating  the  expression  for  c/>i  in  equation  (42)  with 


11 


respect  to  z,  recalling  that  ^  =  1//?||,  and  using  equation  (44)  for  l//fy,  we 
find  that 


dz 
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We  can  obtain  an  expression  for  the  rate  of  change  of  02  with  respect  to  z  in 
a  similar  fashion: 
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d(j>2 
dz 
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w     3  J  ,       UWV3  1 
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(49) 


\r?W>  72^ 

+  0(62a), 

where  0^  =  d(f>y  /dz  must  be  determined  by  using  Maxwell's  equations  to 
determine  the  evolution  of  the  optical  field. 

Note  that  the  electron  energy  (7  in  (47)  depends  on  both  the  phases 
0!,02>  whose  evolution  is  determined  by  equations  (48)  and  (49).  Note  also 
that  the  evolution  of  these  phases  depends  on  the  electron  energy  7. 

If  we  now  compare  the  higher  order  terms  in  the  expressions  for  d^/dz, 
d0i/dz,  and  d02/dz,  we  see  that  unless  £  =  0,  oscillatory  solutions  for  7,^1, 
and  02  will  not  exist.  Since  the  existence  of  oscillations  in  7  implies  a  transfer 
of  energy  from  the  electron  to  the  optical  field,  and  this  is  essential  to  the 
operation  of  the  FEL,  we  find  that  we  must  take  £  =  0  in  equations  (47), 
(48),  and  (49). 

At  this  point  we  have  two  independent  small  parameters:  e  and  ea  =  \i 
(as  we  have  not  yet  specified  a).  From  equation  (25)  e  is  a  measure  of  the 
smallness  of  A^  compared  to  L,  while  \TXI2  measures  the  energy  7  of  the 
electron  (see  (28)).  As  in  [12],  near  the  resonant  state  (7  «  7r),  there  are 
two  distinguished  limits  corresponding  to  a  =  0  and  a  =  2/3.  The  latter 
case  is  known  as  "slow  resonance"  and  is  substantially  more  difficult  to  solve 
than  the  case  when  a  —  1,  and  for  this  reason  we  shall  restrict  our  attention 
to  the  a  =  1  case. 
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Equations  (47),  (48),  and  (49)  form  three  coupled  ordinary  differntial 
equations  for  7,  <^,  and  <f>2.  The  coefficients  in  these  ODEs  are  slowly  varying 
functions  of  z  =  tz.  Though  the  single  mode  case  (see  [12])  can  be  cast 
into  Hamiltonian  form,  it  does  not  seem  possible  to  express  the  two  mode 
equations  in  Hamiltonian  form.  One  must  also  study  how  the  optical  wave 
interacts  with  and  evolves  according  to  the  presence  of  the  electron  beam. 
(We  have  derived  the  equations  for  a  single  electron,  but  generally  in  the 
electron  beam  of  a  FEL  there  are  many  electrons). 

In  order  to  maintain  the  coherence  and  monochromicity  of  the  radiation, 
the  evolution  of  the  optical  wave  should  be  slow.  Thus  the  leading  order  parts 
of  the  signal  parameters  by  and  (f)\  must  either  be  constant  or  depend  on 
the  prescribed  slowly  varying  wiggler  parameters.  These  signal  parameters 
in  the  radiation  field  are  determined  by  the  wave  evolution  process  which  is 
governed  by  Maxwell's  equations.  Though  determination  of  the  asymptotic 
behavior  of  the  signal  parameters  is  relatively  straight  forward  in  the  one 
mode  model  (see  [12]  for  example),  it  still  is  not  quite  clear  how  the  wave 
evolution  occurs  in  the  two  mode  model. 

We  leave  the  further  details  of  this  analysis,  and  the  study  of  the  three 
ODEs  in  (47),  (48),  and  (49)  to  a  future  report. 
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